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Abstract 



o 

We show that in the Klein projective ball model of hyperbolic space, the hyperbolic Voronoi 

diagram is amne and amounts to clip a corresponding power diagram, requiring however al- 
\—A gebraic arithmetic. By considering the lesser-known Beltrami hemisphere model of hyperbolic 

geometry, we overcome the arithmetic limitations of Klein construction. Finally, we characterize 
t— I the bisectors and geodesies in the other Poincare upper half-space, the Poincare ball, and the 

Lorentz hyperboloid models, and discusses on degenerate cases for which the dual hyperbolic 

Delaunay complex is not a triangulation. 

u Key words: Hyperbolic geometry; Voronoi diagram; power diagram; amne diagram; Mobius 

& diagram; Delaunay complex; Delaunay triangulation; hemisphere model; Klein ball model; 

i_i Poincare ball model; upper halfspace model; Lorentz hyperboloid model 

(N 
> 

^ 1 Introduction 

<N 

Given a finite distinct point set V of a space X equipped with a distance function d(-, •), the Voronoi 
diagram PDJ of V tessellates[j]the space into proximity regions called Voronoi cells: 

7-H Vor(p) = {x e X d(p,x) < d(q,x), Vq £ V} ,Vp £ V. (1) 

The ordinary d-dimensional Voronoi diagram is obtained by taking the Euclidean distance 

d{p,q) = \\p — q\\e = \/ (p — Q,P — q)e = \/Yli=i(Pi ~ 1i) 2 (square root of the Euclidean inner 
product of the vector difference) . This Euclidean Voronoi diagram has been extensively studied |23j 
in relation with the Delaunay triangulation. The Delaunay triangulation [T3] DT('P) of a point 
set V is defined as a triangulation such that no point of V falls strictly inside the circumscribing 
spheres of its simplices anchored at V . The Delaunay triangulation is unique for points in general 
position^] (no collinear nor co-spherical degeneracies) . 

* Corresponding author. Fax: (+81) 3-5448-4380. Frank.Nielsen@acm.org, Sony Computer Science Laboratories 
Inc., 3-14-13 Higashi Gotanda 3F, Shinagawa-Ku, Tokyo 141-0022, Japan. 

^Richard . Nock@niartinique . univ-ag. fr CEREGMIA UAG, Campus de Schoelcher, BP 7209. 97275 Schoelcher, 
Martinique, France. 

1 Pairwise disjoint interior tiles fully covering X. 

2 That is, no fc-flat containing k + 2 points and no fc-sphere passing through k + 3 points (for 1 < k < d — 1). 
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In this work, we consider the Voronoi diagram in the d- dimensional real hyperbolic geometry 
space X = M d (K,) of constant negative sectional curvatur^K. We refer the reader to the survey |12j 
for a concise description of the five standard models of hyperbolic geometry and their relationships. 
We call those five models: 

• the Poincare Upper half-space (U), 

• the Poincare ball (P), 

• the Klein ball (K), 

• the Lorentz hyperboloid (L), and 

• the Beltrami hemisphere models (B). 

The above B/K/L/P/U naming and mnemonic of those models does not necessarily reflect the 
historical development of non-euclidean geometry (see [21] for a historical perspective of the first 
150 years of hyperbolic geometry). 

The paper is organized as follows: Section [l . 1| briefly recalls prior work concerning the hyper- 



bolic Voronoi diagram in arbitrary dimension, and section 1.2 motivates the study of hyperbolic 
Voronoi diagrams by presenting some applications ranging from computer graphics to computa- 
tional information geometry. Section [2] describes the hyperbolic Voronoi diagram in the hemisphere 
model and show how to compute it using off-the-shelf algorithms relying only on rational arithmetic. 
Section [3] describes the bisectors and geodesies for the five models, and discusses on degenerate cases 
for which the dual Delaunay complex it not a triangulation. Finally, section [4] concludes this work. 
Appendix [A] reformulates our former work [22] in arbitrary dimension for sake of completeness. 



1.1 Prior work 

Since we consider the general d-dimensional setting, let us recall only prior work tackling the 
hyperbolic Voronoi diagram constructions in arbitrary dimension. Boissonnat and Yvinec [10J (pp. 
449-454) proved that the complexity of the hyperbolic Voronoi diagram of a n- point set V <E H d (-1) 
is 0(n^l ) using the Poincare d-dimensional upper half-space model (U). They proceed by exhibiting 
an injective correspondence between the Euclidean and the hyperbolic diagrams using two successive 
projections. Nielsen and Nock [22] showed that the hyperbolic Voronoi diagram in the Klein ball 
model (K) amounts to compute a clipped power diagram [lj. Appendix [A] recalls this construction, 
extending [22] to arbitrary dimension. Bogdanov et al. [5J proved that the hyperbolic Delaunay 
triangulation (dual graplj^Jof the hyperbolic Voronoi diagram) in the Poincare ball model (P) can 
be obtained by removing the simplices of the Euclidean Delaunay triangulation intersecting the 
bounding sphere. Their approach requires only rational arithmetics^] 

In this paper, we study the hyperbolic Voronoi diagram using the lesser known Beltrami hemi- 
sphere model [12] (B) and hyperboloid model (L), and show that the hyperbolic Voronoi diagram 
in those models amount to compute an affine diagram. 

3 For spherical geometry, the Gaussian curvature of the d-dimensional sphere of radius r is re = \. The curvature 
tends to (flat) when r — > oo. 

4 The Delaunay complex is however not always a triangulation when points are in degenerate positions. See 
Section [3] 

5 However, we need to stick to the Poincare model since conversion formula may introduce square root operations. 
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1.2 Applications 



Hyperbolic Voronoi diagrams find applications in various fields of computer science. In computer 
graphics, Walter [25] proposed an hyperbolic image browser with its user-friendly interactive inter- 
face. In network, Tanuma et al. [24] extend the greedy geometric routing of Kleinberg |19j using 
the dual structure of hyperbolic Voronoi diagrams. In information geometry [18], the Riemannian 
geometry induced by the Fisher information matrix of location-scale families of probability distribu- 
tions (e.g., normal, Cauchy or Laplace distributions) amounts to hyperbolic geometry M d (K), where 
the curvature constant k depends on the density shape of the location-scale family. It follows that 
the statistical Voronoi diagrams with respect to the Fisher-Rao Riemannian distance [18] amounts 
to compute hyperbolic Voronoi diagrams. 



2 Hyperbolic Voronoi diagram in the Beltrami hemisphere model 

The (i-dimensional hyperbolic space can be studied using various analytic models embedded in W. d+1 
using an extra dimension [T2J denoted by xq. The hemisphere model S + C of hyperbolic space 

with constant curvature k = — 4^ is defined on the positive half-sphere of radius r = \l 



x e R d+1 



d \ 

x 2 = r 2 = and x > \ . (2) 



K 

i=0 



Although less prominent in the literature, this hemisphere model was first reported by Bel- 
trami [3] in 1868. The hyperbolic distance between two points p and p' on S + is expressed [20] in 
the hemisphere model by: 

d B (p,p f ) = arccosh (l + /"^'f^ ) , (3) 

i . 



where point a denotes the southern pole of the sphere of radius r = y — -: a = (— r, 0, 0), 

arccosh(x) = log(x + Vx 2 — 1) for x > 1, and (p,p') E = Yli=oPiPi ^ s ^ ne Euclidean inner prod- 
uct. We use the subscript E to distinguish this inner product from the Lorentzian inner product 
introduced in Section [31 

It follows from Eq. [3] that we have: 

co S h(d B (p,p / ))-l= 1 ~ 2 (P, f , (4) 

r z p p' 

where cosh(a;) = eX+ ^ is the hyperbolic cosine function. Thus the bisector Bis(p,p') of two points 
p and p' on the hemisphere model is thus expressed by: 



Bis(p,p') 



1 - (p, x )e _ 1 - ¥, x )i 



r 2 p x 



r 2 p' xo 



i=l 



Pi Pi 



1 \ 



+ 



1 1 



Po PoJ Po Po 



T = 0, i £ 5 



(5) 
(6) 
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Observe that the bisector relaxed to lR. d+1 is a hyperplane with the coefficient corresponding to 
the term xq vanishing (that is, a vertical hyperplane). Therefore the bisector on the hemisphere is 
a vertical (d — l)-sphere. 

To compute the hyperbolic Voronoi diagram on the Beltrami hemisphere, we can thus choose 
hyperplane Hq : xq = 0, compute the affine diagram, and clip it with the ball Ya=i x f < f 2 - We 
then either lift vertically this diagram onto the hemisphere, or manipulate the diagram internally 
using its affine representation on Hq. (For example, point location can be handled directly on Hq.) 

On hyperplane Ho, we write the bisector equation as: 

Bi B (p,p'):(xX--) +-~\ = 0, (7) 
\ Pq Po/ E d Po P 

where (•,•)#<* denotes the d-dimensional Euclidean dot product on coordinates x\,...,x^. Affine 
diagrams can be built equivalently using power diagrams [lj: A power diagram for a given a set 
B = {Bi, B n } of n balls (with Bi = Bahg^Cj, Wi = rf) for i G {l,...,n}), is defined as the 
minimization diagram [§] of the corresponding n functions: 

Di(x) = (a - x,a - x) E - Wi, (8) 

where Wi is the weight associated to center point q. The power bisector Bn{Bi,Bj) of any two 
distinct balls Bi = Balled, Wj) and Bj = Ball£;(cj, Wj) is the radical hyperplane [1] of equation: 

Bi n (Bi, Bj) : 2(x, Cj - q) e + (q, Cj) E - (cj, cj) E + wj - Wi = 0. (9) 

The last equation shows that power diagrams are additive weighted (ordinary) Voronoi dia- 
grams |10| . 

To find the equivalence of the Beltrami hemisphere affine diagram with a power diagram, we 
map the points p to equivalent balls B with centers c and weights w such that: 

p 1 , / \ 1 , (p>p)e 

C=— , W= h(c,c) E = h 2 ■ 

2po Po Po <±Po 

This calculation requires only rational arithmetic but input shall be given using d + 1 coordi- 
nates. (Klein hyperbolic Voronoi diagram used (^-coordinates but requires algebraic coordinates, 
see Appendix [A]) 

Although the worst-case combinatorial complexity 0(nJ^) was already obtained in the Klein 
ball model [22], the Poincare ball model [5], and the Poincare upper space model [10], this 
novel hemisphere construction allows us to use off-the-shelf power diagram constructions to build 
the hyperbolic Voronoi diagram using only rational arithmetic. Note that since the hyperbolic 
Voronoi diagram is a particular case of the Riemannian Voronoi diagram, it follows that the worst- 
case infinitesimally-scaled Euclidean Voronoi diagram (from equivalent polytopes on the d + 1- 
dimensional curve of moments [TU]) can be obtained. Therefore the complexity of the hyperbolic 
Voronoi diagram is O(n^l). 

3 Voronoi bisectors and geodesies in the five standard models 

Hyperbolic geometry can be studied under the framework of Riemannian geometry [11] . Table [l] 
summarizes the Riemannian metrics and the Riemannian distances for the Beltrami hemisphere 
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Model 


Riemannian metric ds 2 


cosh^=, ci(-,-) : Riemannian distance 


Klein (K) 


( ds E i i x >dx) E \ 


1-<P,9>b 


K \ i UtII 2 ' n imi 2 ~\2 1 
\ i — IWIe \ i -~\\ x \\e> J 


//■-I ll i? \ /i ii ii9n 

V(i-lbll|)(i-lkll|) 


Poincare (P) 


4ds| 

My / -i II no \o 

(i-IMI|) 2 


-I . 2||p- g ||| 

-L H nO\/-i ii nO\ 
(HMIrOCHMIb) 


Upper (U) 


— K — 

x d 


i | lb-9ll| 


Hyperboloid (L) 


fcdx 2 , + ds|. 


<P.9>L 


Hemisphere (B) 




1 | 1-(P.9>E 
Kp<)<?0 



Table 1: Models of real hyperbolic geometry M d (n) of constant section negative curvature k as 
Riemannian geometries: Riemannian metrics and distances. All models are conformal (angle- 
preserving) except Klein model that is not conformal (except at the origin). 



(B), the Poincare upper half-space (U), the Poincare ball (P), the Klein ball (K), and the Lorentz 
hyperboloid (L) models. The five model domains are (r — - - 



V 



{(xo,r,x 2 , —Xd) | x > 0} (upper space), 
< (0,xi,...,x d ) 



K 



(x ,xi, ...,Xd) 



d 1 

x 2 < r 2 > (Poincare ball), 
i=l J 

d \ 

xf < r 2 and xo > > (Beltrami hemisphere), 

i=0 J 



(x ,xi 



E 

i=i 

E 

i=i 



(Klein ball), 



(10) 
(11) 

(12) 
(13) 



2 _ 2 
Xj x 



-r 2 and xq > > (Lorentz hyperboloid). (14) 



The Poincare and Klein ball models and the Poincare upper half-space model have been often 
considered from an application point of view since their domain dimension does not need to increase 
the ambient dimension and they can be easily displayed on a computer screen. 

These five models are all linked in R d+1 as described in [12] (Figure 5, page 70). The maps from 
the hemisphere to the other models are either central or vertical projections explicited in [12]. Let 
us summarize the bisector expressions in those five standard models: 



BiflO, q) 

Bii<-(p, q) 
Bi P (p, q) 



E< 

i=i 
x, 



* _ Ei 
qo po 



l 



+ 



l 



l 



Po % 



p¥e-p\Ji - UVe j + V i " U\e ~ V 1 " Me = °> 

(x,x) E I - 11^112 \ ||„||2 ) + 2( X, - ; ; — — : r; n-r ) + 



We 1-klle 



1 



|2 



i- HI/ e 



(15) 
(16) 
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Model 


Geodesic 


Bisector 


Hemisphere 


arc of circles 


vertical spherical portion of the 
hemisphere 


Klein ball 


line segment 


hyperplane 


Poincare ball 


circular arc perpendicular to the 
bounding sphere, or straight 
lines passing through the origin 


spherical portions, or hyper- 
planes passing through the origin 
when ||pl| = \\q\\ 


Upper half-space 


arc of circles or vertical lines 


spherical portions or hyper- 
planes 


Hyperboloid 


hyperbola 


intersection of hyperplane pass- 
ing through the origin with the 
hyperboloid 



Table 2: Hyperboloic Voronoi diagrams: Characterization of geodesies and bisectors for the five 
standard models. 



II Iklll 



i-IWI 2 i-NI 2 " ' (17) 

1 \\P\\E 1 \\Q\\E 

d 

Bi L (p,q) : (x, q - p) L = (p - q )x + - Pi)x, = 0, (18) 

i=i 

1 1 A , „/ Q P \ . (p>p)e (Q,Q)e 



Biu(p,q) : {x,x) E [ + 2( x, ^ - ^ ) + ^™ _ v*l^ = . (19) 

\Pd Qd) \ Qd Pd/ E Pd Qd 

(20) 

where (x, y) h = — so2/o+Sf=i x iVi denotes the Lorentziarj^] inner product. Observe that the Lorentz 
bisectors relaxed to M. d+1 are hyperplanes all passing through the origin. (The geodesies are hyper- 
bola arcs contained in a hyperplane passing through the origin.) The Lorentz hyperbolic Voronoi 
diagram can thus be obtained as the intersection of a {d + l)-dimensional affine diagram with the 
hyperboloid upper sheet L + . For odd dimension d = 2k + 1, we can thus build optimally the 
hyperbolic Voronoi diagram in the Lorentz model since O(nW) = O(ralYI) = 0{n k+1 ). That 
is, the extra dimension xq for odd parity does not exhibit the dimension gap in the combinatorial 
complexity of the Voronoi diagram. The hyperboloid model has been used by Galperin |16] to pro- 
pose the model centroid in a closed-form expression. (Indeed the hyperboloic Karcher centroid |17| 
does not admit an analytic expression.) 

Table [2] characterizes the nature of the bisectors (bisecting sites) and geodesies (linking sites) 
for those models. 

It is well-known that for points in general position (no collinear nor co-spherical points), the 
dual of an ordinary Voronoi diagram is the unique Delaunay triangulation [10] . This duality holds 
for the hyperbolic Voronoi diagram and the hyperbolic Delaunay triangulation (with the empty 
sphere property), provided general position of the sites. Hyperbolic Voronoi diagrams can be fairly 
different from ordinary Voronoi diagram. For example, Figure [2] gives an example of the hyperbolic 
Voronoi diagram when all but one site are co-circular and close to the bounding circle with a single 
site centered at the ball center. In that case the dual Delaunay complex is not a triangulation 
although unique (a (n — l)-ary tree of depth 1). We refer to [7] for a stability analysis of the 



3 Although attributed to Hendrik Lorentz, the model is due to Karl Weierstrass [15] . 
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Figure [T] depicts some degenerate cases in the ball models for co-spherical points centered at the 
origin, and their corresponding diagrams on the upper half-space model. 

Klein ball Poincare ball Upper halfspace 




Figure 1: Degenerate hyperbolic Voronoi diagrams: (a): In the Poincare/Klein ball models when 
point radii are all equal (wheel diagram with bisectors passing through the origin.) (b): In the 
upper halfspace when the d-th. coordinate are all equal. 
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Figure 2: A hyperbolic Voronoi diagram with all cells unbounded (n = 9). The dual graph 
corresponds to a tree: (a) Klein, (b) equivalent power diagram, (c) Poincare ball, and (d) Poincare 
upper halfspace. 

Delaunay triangulation on generic manifolds. 

4 Conclusion 

In the Klein ball model, we extended the planar case j22j and proved that the hyperbolic Voronoi 
diagram of a finite point set in arbitrary dimension is affine. It followed a simple method to 
optimally calculate the Klein hyperbolic Voronoi diagram by computing a corresponding power 
diagram clipped to the bounding sphere. This method however required non-rational arithmetic 
(square root operations in the denominator) on <i-dimensional coordinates as noticed in [5]. To 
overcome this drawback, we described a novel approach in the hemisphere model, and showed that 
the hyperbolic Voronoi diagram also amounts to compute an affine diagram using only rational 
arithmetic on (d + l)-dimensional coordinates. From the view point of computational geometry, 
those various analytical models of hyperbolic space allow one to take advantage of the merits of each 
model. The Poincare/Klein ball models are both conformaj^] at the origin with geodesies passing 
through the origin being straight lines. For geodesic walking-type algorithms like computing the 
hyperbolic centroid, Barbaresco [2j considered the Poincare ball model and perform a hyperbolic 
translation at each step to set the current point to the origin in order to ensure that geodesies 

7 Preserving angles. That is, the metric is locally a scaled Euclidean metric (identity matrix). 
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departing from that current point are straight line segments. This could have been also performed 
using the Klein model. When manipulating these models, we may perform a hyperbolic rigid 
motion [20] to choose a convenient origin [4jj^] 
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A Hyperbolic Voronoi diagrams in the Klein projective ball model 

This appendix extends |22j to arbitrary dimension. In the Klein ball model /C(k), the hyperbolic 
geometry of curvature k = — ^ is embedded [11] inside a Euclidean ball IB = {x E M rf | \\x\\e < r} 
of radius r with the corresponding Klein hyperbolic distance dx(p,q) between two points p and q 
expressed by: 

d K (p, q) = x arccosh ( = = == | , (21) 

V^(l-b|||)(l- 

where arccosh(x) = log(x + \/x 2 — 1) for x > 1. 

Consider the Klein bisector Bix(p, q) of any two distinct points p and q: 

^k(p,q) = {x G JC(k)\ d K (p,x) = d K (q,x)} , (22) 

Since arccosh(x) is a monotonic function preserving the distance order and \J— k is a multi- 
plicative constant, it follows that: 




Bi K (p,q) = < x G 



l- WD(i-INH) J(i-\\ q f E )(i-\\ x \\ E y 



(23) 



\E/\ ~ H^llS/ V \ x Wi\\E 
Therefore, the bisector Bx(p,q) is the loci of the points x satisfying: 



^V 1 - ll* -py 1 ~ ^ ) + v 1 - ^ - v 1 - = o- (24) 

This is a linear equation (a, x) E + b = in cc, with: 

l-\\pf E - Py /l-\\ q \\%, 6=^l-|k|||-^l- W|. (25) 

That is, Klein bisectors are hyperplanes. It follows that the Klein hyperbolic Voronoi diagram 
is an affine diagram [9] with all its Voronoi cells convex Affine dia grams can be universally built 
from an equivalent power diagram [9]. The power distance dn(B,x) of a point x to a Euclidean 
ball B = Balls (c, r) of center c and radius r is defined as: 

du(B, x) = (c — x, c — x) E — r 2 . (26) 

Given a set B = {Bi, B n } of n balls with Bi = Balled, n), the power diagram is defined as 
the minimization diagram [9] of the corresponding n functions: 

Di(x) = (ci - x,Ci- x) E - rf . (27) 

The power bisector B]j(Bi, Bj) of any two distinct balls Bi = Balls (cijl) an d Bj = Balls (cj,i"j) 
is the radical hyperplane [1] of equation: 

Bin(-B», Bj) : 2(x, cj - q) e + (a, cf) E - {cj,Cj) E + r) - rf = 0. (28) 



9 It is a particular affine diagram since all its cells are guaranteed non-empty. 
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This power bisector linear equation shows that power diagrams are also affine diagrams. Let 
us report the mapping that associates to points in the Klein ball corresponding Euclidean balls 
so that their bisectors coincide pairwise. By identifying the bisector equation Eq. [24] with that of 
Eq. [9J we find that point p% is mapped to a ball B{ of center a with squared radius rf, where: 

d = c(B i ) = — F M , rhr 2 (B^A_ 1 (29) 



2,/1-INIl' ' ' 4 1 " lift III) Ji-WnWl 



Observe that the ball centers may span the full <i-dimensional Euclidean space Mr. The radii may be 
negative and are interpreted as balls with imaginary radii (or as weighted points with real weights, 



see Eq. 26). Solving for the quadratic expressiorj | of Eq. 29, we find that r (pi) < if and only if: 

\\pi Hi < 4(V5 - 2) ~ 0.944272. (30) 

Thus to compute the Klein hyperbolic Voronoi diagram of point set V = {p\, ...,p n }, we first 
construct a set of corresponding balls 

B = {B 1 = Ball E (c(pi),r(pi),..,S n = Ball B (c(p n ), r(p n )}, (31) 

compute the power diagramp] of £>, and consider the restriction of this diagram to the unit ball B d . 
Figure [3^a) and Figure [3^b) depict the Klein Voronoi diagram and corresponding power diagram, 
respectively. It follows another proof of the combinatorial complexity of the hyperbolic Voronoi 
diagram using the Klein model instead of the Poincare upper halfspace [10j (pp. 449-454) or the 
Poincare ball model [5]: 

Theorem 1 The hyperbolic Voronoi diagram of n d- dimensional points can be obtained in the 
Klein ball model as an equivalent power diagram clipped with the interior of the unit ball. The 
diagram has combinatorial complexity 0(n^^) and can be built in optimal 0(n log n + n^l) time. 

Figure [3](a) displays the Klein hyperbolic Voronoi diagram built from an equivalent power 
diagram (Figure [3^b)). Note that all bounded Voronoi cells necessarily fall inside the ball, and 
therefore the construction can be made output-sensitive using Chan's algorithm |13j . 

The following section describes the mapping functions to convert from the Klein model et the 
Poincare ball, upper halfspace and Lorentz models. 



B Converting Klein model to other models 

We report formula 12 and their inverse for converting from/to Klein coordinates into the Poincare 



ball (P), upper halfspace (U) and Lorentz coordinates(L): 



10 Let I — \\p\\% £ [0, 1]. We seek when 4 ( 1 '_ i ) — 775=1 < 0. That is, we solve for the root of equation / — 4Vl — I in 
[0, 1]. This is a quadratic expression I 2 + 161 — 16 = which admits a unique root in [0, 1]: 4(\/5 — 2). 

11 Power diagrams can be computed in arbitrary dimension as convex hulls or halfspace intersections (unbounded 
polytope) using free software like QHull (http://qhull.org), CGAL (http://www.cgal.org) or LEDA (http://www. 
algorithm ic- solutions . com/leda/[ |. 

2 Wlog., we assume k = — 1. Otherwise, we first rescale by \J 



12 



Figure 3: (a) The Klein Voronoi diagram calculated from the restriction of the power diagram (b) 
to the unit ball. Converting the Klein diagram into the Poincare Voronoi diagram (c) and the 
upper half-space Voronoi diagram. Observe that the combinatorial complexity of the cells does not 
change (n = 16 points). The power diagram may have cells empty of sites and ball centers may be 

located outside the square of side length r = J — K 
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